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Abstract Some charmonium-like resonances such as X(3872) can be interpreted as possible D'*'D^ molec- 
ular states. Within the quark model, we study the structure of such molecular states and the similar 
molecular states by taking into account the light meson exchange (ir, 77, p, uj and a) between two light quarks 
from different mesons. 
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1 Introduction 



Some recently discovered narrow charmonium- 
like resonances have aroused great theoretical in- 
terest. The typical X(3872) state, which was dis- 
covered by the Belle Collaboration |l| and subse- 
quently confirmed by the CDF Collaboration Q and 
BABAR Collaboration [3|, is summarized M x = 
3872.2 ±0.8MeV, T x = 3.0±^ ± 0.9MeV Previ- 
ously its quantum numbers were inferred most proba- 
bly J PC — 1 ++ , while the masses of the corresponding 
charmonium states in the quark model 2 3 P 1 or 3 3 P 1 
are hundreds MeV above M x . Recently, a study of 
the ir + Tr~ir° mass distribution from the A(3872) de- 
cay by the BABAR Collaboration favors a negative 
parity assignment 2 _+ 5]. However, still the mass of 
the charmonium candidate in the quark model 1 X D 2 
disagrees with M x . 

Since it is difficult to assign X(3872) to any con- 
ventional cc charmonium state in the quark model, 
other alternative explanations prevail. Because M x is 
close to the D*D threshold, X(3872) was interpreted 
as a loosely bound molecular D*D state [B-S] or such 

a molecular state with some admixtures of pJ lib and 

r— 1 ' 
uiJ/ip uM- It was also suggested as a tetra quark 

state dominated by a diquark-antidiquark structure 

11 1 or a ccg hybrid state The proximity of the 

X(3872) to D*D threshold may also imply that the 

cusp scenario should be treated seriously [13f . 



The D M D M molecular states were proposed 
To bind the two mesons 



14, 



15 



many years ago 
together in a molecular state, additional interaction 
beyond the quark potential model should be intro- 
duced between the two mesons. Naturally the me- 
son exchanges were first considered. The inter-meson 
potentials from light meson exchanges can be eas- 
ily treated in the frame of the heavy quark effecitive 
theory (HQET). The studies showed that the cut- 
off parameters in the form factors are critical to the 
formation of the molecular states Other 
interactions such as the gluon exchanges were also 
considered [13, [ll, ■ 

Basically, the form factor is related to the in- 
ner quark wave functions of the meson state. To 
avoid the uncertainty from the form factor, in this 
work we will investigate the possible molecular states 
directly using the quark wave function within the 
quark model. Our calculation is based on the multi- 
Gaussian function expansion of the quark wave func- 
tion of the molecular state, which is a rather simple 
and efficient variational method to study many-body 
ground bound states [2CH23|. 

We will assume that the light mesons (7r, 77, p, w 
and a) are exchanged between the light quarks from 
different hadrons. The mechanism can be understood 
from the chir al q uark model 24j which was proposed 
by Weinberg [25[ and formulated by Manohar and 
Georgi 24 1. The key feature of the chiral quark model 
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is that between the QCD confining scale and the chi- 
ral symmetry breaking scale, QCD can be roughly 
described by an effective Lagrangian of quarks and 
pseudo-scalar mesons of Numbu-Goldstone bosons. 
The strong interaction between hadrons is thus de- 
scribed mainly by the exchange of pseudo-scalar 
mesons at long-range. In practice, the chiral quark 
model can be further extended with more mesons to 
account for the intermidiate-range interaction. The 
light vector resonances p, K*, etc. can be realized as 
the dynamical gauge bosons of a hidden local sym- 
metry [26| . The scalar meson a can be explicitly in- 
troduced from the linear realization, as in the linear 
sigma model. Chiral quark model was widely used in 
the study of nuclear forces 27, 28j and recently used 
to study the molecular states [29J. 

In section 2, we will introduce the quark model 
with the meson exchange interactions in our work. 
In section 3, the multi-Gaussian function expansion 
method and the configuration space of the molecular 
states are presented. We will study both the 
and molecular states. Finally we will give a 

short summary. 

2 The quark model with meson ex- 
change 

The Hamiltonian in a quark model can be writ- 
ten as 

h = z^+ly-i e ■ w + f > • p ^ ■ s ^ ■ 



i<3 



(1) 

The first part is the non-relativistic kinetic energy, 
where m/s are the constituent quark masses. The 
second part is the central potential, where Fi = 
|Ai are the well known SU(3) C Gell-Mann matrices. 
Apart from a constant, the central potential usually 
is a combination of the one gluon exchange coulomb 
potential and the linear confinment: 



V c = -- 



■Mo, 



(2) 



where r\j = |t\ — r 3 -| is the distance between quark 
i and quark j. The last part is the color-magnetic 
interaction, where S , i =|<Tj are the Pauli matrices. 



In our work we will use the Bhaduri quark model 
which is a rather simple non-relativistic quark poten- 
tial model [30]. The potential function of color- 
magnetic interaction reads 



V s = — l — e - r »' r ° 

m x rrij rgry 
The model parameters are 

k = 102.67MeV, a = 0.0326(fm/McV) 2 . 

M = 913.5MeV, r = 0.4545fm, 

m u =m d = 337MeV, m s = 600MeV, 

m c = 1870MeV, m b = 5290Mev. 



(3) 



(4) 



The quark model is very successful in describing 
the hadron properties at low energy under the QCD 
confinemnet scale. As is well known, due to the color 
symmetry, the above quark potential does not pro- 
vide direct interaction between hadrons. Here we will 
consider the exchange of the light mesons between the 
two light quarks from different hadrons based on the 
chiral quark model. At present, we will consider the 
exchanges of 7T, rj, lj, p, and a mesons only. The 
interaction Lagrangian densities are 



1. Pseudoscalar 



2. Scalar 



C a = g a ip(x)?p{x)(f){x), 



(5) 



(6) 



3. Vector 



C v = g v ip{x)-f ti il)(x)v ti (x) 

fv 



2m, 



(7) 



Here m q is the constituent quark mass, ip(x) is the 
constituent quark Dirac spinor field. f{x), </>(x), 
v' 1 (x) are the pseudoscalar, scalar and vector interme- 
diate meson fields respectively. In the case of isovec- 
tor mesons, the above meson fields should be replaced 
by T-<p, T-(p or t-u m respectively, where r are the iso- 
spin Pauli matrices. From the effective Lagrangian we 
can obtain the effective potential between two quarks. 
In the coordinate space |3lj |. 



I. Pseudoscalar 

V P (r) = j^™ P (^YrM^ .«x 2 , (8) 
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2. Scalar 



V.(r)-. 



47T 



Y(m a r)- 



3. Vector 



.11 f m A 2l V f \ 

+ —m v [ — -Y 1 (m v r)cr 1 -cr 2 . 
47r \m q J 6 



m„ \ 1 



-Yi{m s r) 



9vf v ( rn % 
— — m v — 
47r \ m c 



1 , . 1 . 

- Yi {m v r) + - Yi (m v r)(T 1 ■ er 2 



(9) 



(10) 



We have neglected the momentum dependence of 
all potentials in the present work. The spin-orbit in- 
teraction and the tensor-force are also dropped out 
here since we only consider the molecular ground 
states. The functions Y, Y 1 are defined as follows: 



Y(X): 



YAx)=- 47n5 3 (a;). 

x 



(11) 
(12) 



Note that on the quark level we need not introduce 
the form factors to treat the ^-interaction. The me- 
son's quark wave functions will naturally smear out 
the singularity. 

The qq potentials are obtained from the above qq 
potentials by a G-parity transformation. Also there 
should be an isospin factor t x -t 2 in the above effective 
potential if an isovector meson is exchanged. 

From PDG0, the masses of exchanged mesons 
are taken to be 



= 134.98MeV, 
772 p = 775.5MeV, 
777 CT = 600.0MeV. 



m„ = 547.85MeV, 
m,, = 782.7MeV, 



The masses of relevant heavy flavor mesons are: 



m D o = 1864.5MeV, 
m Dt0 = 2006.7MeV, 
777 S o = 5279.5MeV, 
777 B , = 5325.1MeV. 



777 D ± = 1869.3MeV, 
m D .± = 2010.0MeV, 
to s ± =5279.1MeV, 



The coupling constants are taken from Ref. Q: 

g„ = 3.30, 



9, = 2.74 
.9p = 3.46, 
9„ = 5.28, 



9n = 2.05, 
/p = 1-45, 

L = o, 



which are extracted from the mcson-nucleon coupling 
constants in the well-known Bonn mode l 13 ll using 
the single-quark operator approximation [32j. 



3 The heavy flavor molecular states 

In this work, we consider the possible molecu- 
lar states constructed from the pseudo-scalar heavy 
mesons (D, B) and their vector partners (D* , B*). 
The states involving strange mesons are not consid- 
ered here. The corresponding charmed combinations 
are: D-D, D*-D* , D*-D. Since the charmed mesons 
belong to the I = | representation of isospin SU(2), 
the possible isospins of the D^-D^ system are 1 = 0, 
1. Following Ref. 0, we label the D-D, D'-D and 
D*-D* systems as <&* u and $j} respectively, 

while the D*-D system with negative charge conju- 
gate (C = — 1) is labeled as $ 7J . Below we pick up 
the neutral state to represent the isospin multiplet: 



1. 1 = 



& 00 = -^(D°D0-D-D+) 



(13a) 



%? = — (D*°D*° - D*-D*+) T (J = 0, 1, 2) 

(13b) 

(13c) 



^01 



-.±[D°D*°-D + D*- 
2 L 

.C(D°D*°-D-D*+)], 



2. 1=1 

*!„ = 



(D°D° + D-D^ 



1 

1 + D *- D *+\ 

v2 



= ~[D°D*° + D+D— 
-C{D°D* a + D-D*+)]. 



(13d) 

(J = 0,l,2) 
(13e) 

(13f) 



The states of the B^ — B^ combinations are con- 
structed similarly but named after fi. 

To calculate the molecular state, we use the 
Raylcigh-Ritz variation principle. The test wave 
function is taken to be a series of Gaussian functions 
with various widths 21- 23j. In our case of the QqqQ 
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molecular state, the test wave function between the 
Qq and qQ mesons is expanded to [2o| 



^1234 (f 12 ! r 34 : ^1234 ) 



^2a' 1234 ip 12 (r 12 )ip 34 (r 34 ) 

i 

xexp(-/3" 2 3 4 r? 234 ), (14) 



where r u r 2 , r 3 and r 4 are the coordinates of Q, q, q 
and Q, respectively, — Ti—Tj. ; 



is the distance 



between the two meson clusters 

_m Q r 1 + m q r 2 m q r 3 +m Q r 4 



m Q +m q 



m q +m Q 



(15) 



Each of the meson wave functions iptj (r^ ) is also 
taken to be a Gaussian function series 



^(ry) = ^<exp(-/3*r^)- 



(16) 



Numerically the above wave function is deter- 
mined by the variational method in two steps. First 



the meson wave function (| 10|) is determined from the 
potential quark model. Then the wave function (|14[) 
of the molecular state is obtained from the meson ex- 
change potentials with the meson wave functions 
fixed. 

To reduce the amount of computation, the pa- 
rameters ft 1 and a 8 in each Gaussian function series 
are determined also in two steps by one-dimensional 
minimization. We first search a central (5 value using 
a single Guassian function. Then a set {/?'} of 27V+1 
elements is generated by scaling the /3 value up and 
down by a scale factor s [2II ]: 



(17) 



where i = 0,l,...,2iV. 

The bound energies and the mean squared radii 
(rms) (7*1234) 1,72 of the £)(*)£)(*) molecular states are 
listed in Table [U 



Table 1. The bound energies of the D^*'D^ molecular states. In calculation II, all the meson coupling 
constants are reduced by a factor of 0.7, except the tt meson. 



T = 


*00 


*Si 




*Si* 


*S2 


rms(fm) 


I 

II 


1.34 
2.60 


1.37 
2.47 


0.94 


1.28 
3.14 


1.48 
2.54 


£(MeV) 


I 

II 


-29.7 
-1.9 


-30.7 
-2.6 


-45.1 


-24.6 
-1.2 


-24.6 
-2.3 


T=l 




*ii 


*10 






rms(fm) 


I 

II 


1.24 


1.12 


1.04 


1.19 


1.71 


_B(MeV) 


I 

II 


-13.7 


-18.3 


-22.6 


-15.6 


-5.7 



The calculation shows that the meson exchange 
interaction is strong enough to bind the molecular 
states. The typical $o X is the candidate molecular 
state for the X(3872) . However the bound energy of 
$qj is 30.7 MeV, which is too larger than what we ex- 
pect, for X(3872) should be a loose bound molecular 
state. 

However, there is some uncertainty in the meson 
coupling constants on the quark level. In the Bonn 
model for neucleon interaction, the form factors 31 1 

'Al-m 



F a (k 2 )- 



(18) 



are also introduced in the description of the meson 
baryon couplings. Clearly from the meson mass de- 



pendence in the form factor, the effective meson cou- 
pling constants decrease as the mass of intermediate 
mesons increase. 

Next we try to decrease the coupling constants 
to the mesons er, p, uj and 77 with heavier masses by a 
factor A following Ref [9(. The numerical results for 
the case of A = 0.7 are shown in Table [1] as calcula- 
tion II. Now the bound energy of is only 2.6 MeV 
and the rms is 2.47 fm, which meet the interpretation 
of X(3872) as a loose bound molecular state. Other 
possible molecular states left are iso-scalar $ 00 (0 ++ ), 

<I>;*(1 ++ ) and <&; 2 *(2++). 

The similar calculation to B^B^ molecular 
states are shown in Table [5] 



Table 2. The bound energies of B^B^ molecular states. In calculation II, all the meson coupling constants 
except the 7r meson are scaled by a factor of 0.7. 
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T = 


f^oo 


U 01 


"oo 


U 01 


^02 


rms(fm) 


I 

II 


1.13 
1.39 


1.16 
1.32 


0.65 
1.11 


1.00 
1.38 


1.21 
1.42 


£(MeV) 


I 

II 


-46.0 
-13.1 


-49.2 
-17.3 


-81.0 
-11.7 


-44.6 
-9.5 


-44.0 
-15.3 


T = l 


fiio 




"10 


on 


n 12 


rms(fm) 


I 

II 


0.72 
1.01 


0.66 
0.90 


0.62 
0.84 


0.68 
0.93 


0.84 
1.29 


_B(MeV) 


I 

II 


-38.6 
-10.4 


-48.9 
-14.8 


-59.3 
-19.3 


-46.6 
-13.7 


-25.6 
-5.2 



As we expect, the bound energies of _£?(*) 
molecular states become larger. Even when the cou- 
pling constants are weakened by the scale factor 
A = 0.7 in calculation II, the meson exchange is still 
strong enough to bind the molecular states for all 
combinations. 

4 Summary 

Based on the meson exchange interaction be- 
tween light quarks, we have investigated the heavy 
molecular states in quark model. The molecular 
states are described by the four quarks wave function 
expanded as a series of Gaussian functions. The nu- 
merical results show that the light meson exchanges 
of 7r, rj 1 p, uj and a between the light u, d quarks are 
strong enough to bind the heavy molecular states. 

However, the bound energies of the molecular 
states are tens of MeV (up to 80MeV in 
cases) if we adopt the meson-quark coupling con- 
stants from the meson-nucleon coupling constants 
simply using the single-quark operator approxima- 



tion. The results are unreliable as the bound energies 
are somehow close to Aq CD while only the long-range 
meson exchanges are considered in our calculation. 

After we consider the uncertainty of the coupling 
constants which are deduced from the Bonn potential 
of neucleon interaction by decreasing the rj, p, uj and 
a couplings by a factor of 0.7, the A(3872) is well in- 
terpreted as a loose molecular 1 ++ state. The calcu- 
lation also shows that other possible molecular states 
such as 1 ++ and 2 ++ may exist on the threshold of 
D*D* (The ++ scalar state on the threshold of DD 
is complicated from the scalar admixture). Since the 
bind energies of B^B^W molecular states increase 
with the the mass increase of heavy favor, there are 
more such molecular states near the B^B^ thresh- 
olds as we expect. 

The main uncertainty in the work is on the es- 
timation of meson-quark coupling constants. To ob- 
tain a set of reliable coupling constants, we can use 
the model to study the nucleon-nucleon interaction in 
the future. 

We would like to thank professor Shi-Lin Zhu for 
useful discussions. 
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